Coexistence of Weak and Strong Wave Turbulence in a Swell Propagation. 
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By performing two parallel numerical experiments — solving the dynamical Hamiltonian equations 
and solving the Hasselmann kinetic equation — we examined the applicability of the theory of weak 
turbulence to the description of the time evolution of an ensemble of free surface waves (a swell) on 
deep water. We observed qualitative coincidence of the results. 

To achieve quantitative coincidence, we augmented the kinetic equation by an empirical dissipation 
term modelling the strongly nonlinear process of white-capping. Fitting the two experiments, we 
determined the dissipation function due to wave breaking and found that it depends very sharply on 
the parameter of nonlinearity (the surface steepness). The onset of white-capping can be compared 
to a second-order phase transition. This result corroborates with experimental observations by 
Banner, Babanin, Young fl2| . 



INTRODUCTION. 

Wave turbulence is realized in plasmas, liquid helium, 
magnetohydrodynamics, nonlinear optics, etc. A perfect 
example of wave turbulence is a wind-driven sea. The 
major conceptual difference between wave turbulence 
and "classical" turbulence in an incompressible fluid is 
the presence of a characteristic dimensionless parame- 
ter /i, characterizing the level of nonlinearity. Turbu- 
lence is considered to be "weak" if /i <C 1 , otherwise it is 
"strong". In classical hydrodynamic turbulence (i — > oo. 

A more rigorous definition of weak turbulence is the 
following: this is the turbulence which is well described 
by the kinetic equation for waves. These equations are 
the quantum kinetic equations for bosons in the limit 
of very high occupation numbers. They were derived in 
statistical physics in the late twenties [3, [3] and redis- 
covered in nonlinear wave dynamics in the sixties. The 
kinetic equation, describing four-wave resonant interac- 
tion of gravity waves, was named after K. Hasselmann, 
who derived it in 1962-1963 0]. 

The theory of weak turbulence is well-developed 
The kinetic equation has rich families of Kolmogorov- 
Zakharov (KZ) and self-similar solutions, which can be 
efficiently used for explaining a wide range of experimen- 
tal data d, However, today we have a clear under- 
standing of the following fact: even for small values of 
/j,, the theory of weak turbulence may be incomplete. In 
many important physical situations weak and strong tur- 
bulence coexist. 

Even if the weak turbulent resonant interaction effects 
dominate in the greater part of space, strongly nonlin- 
ear effects could appear as rare localized coherent events. 



If they are smooth and regular, they are solitons, qua- 
sisolitons or vortices. However, they could be catas- 
trophic, in which case they are wave collapses, similar 
to self-focusing in nonlinear optics or Lagmuir collapses 
in plasma. Even rare sporadic collapse events can essen- 
tially affect the physical picture of wave turbulence. 

There are two main types of wave collapse events in 
a wind-driven sea. The first is the formation of freak 
waves; this is not a subject of our study. The second, 
which is much more common, is wave-breaking or white- 
capping, which is an essential mechanism of wave energy 
dissipation. It would be hopeless to develop an efficient 
operational model of wave forecasting without an under- 
standing and a proper parametrization of this fundamen- 
tal effect. Meanwhile, a reliable analytical theory of this 
phenomenon is still not developed, while field and labo- 
ratory experimental data are scarce. The most promising 
approach to resolving this problem is a massive numerical 
experiment. 

The most informative experiment would be one that 
could provide a direct numerical solution of the primi- 
tive dynamic equations describing the wave ensemble. In 
1992, Dyachenko, Pushkarev, Newell and Zakharov nu- 
merically solved 2-D focusing NLSE and observed the 
coexistence of self-focusing collapses with weak turbu- 
lence [7!]. Later on, the 1-D MMT (Maida, McLaughlin 
and Tabak) model and its generalizations were solved 
numerically by different authors (see summary in 0). 
Again, the coexistence of wave collapses and weak tur- 
bulence was verified. In our article, we present (as we 
hope, for the first time in the literature) the results of a 
far more detailed experiment. We performed the numer- 
ical simulation of the evolution of an ocean swell using 
two different approaches. 
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In the first, we solved the Euler equations for the 3-D 
potential flow of an ideal incompressible fluid with a free 
surface in the presence of gravity. We used the Hamilto- 
nian form of these equations [loj ]. For gravity waves, 
the parameter of nonlinearity is the average steepness /x. 
We expanded the Hamiltonian in powers of \x up to order 
/i 4 . In the second experiment, we solved the Hasselman 
kinetic equation. 

The comparison of the results demonstrates qualita- 
tive accordance. Both experiments describe expected ef- 
fects, such as the downshift of the spectrum peak, the 
angular spreading of the spectrum and the formation of 
Zakharov-Filonenko spectral tails ~ lo~ a 0, EH- To 
obtain quantitative coincidence of the results, we have 
to augment the Hasselmann equation by an empirical 
dissipation term Sdiss, modelling white-capping effects. 
We tried several versions of this term. The versions 
of Sdiss used in the industrial wave-predicting models 
WAM3 and WAM4 essentially overestimate the dissipa- 
tion for a moderate steepness. The comparison with 
dynamical computations shows that white-capping dissi- 
pation decreases dramatically with decreasing steepness 
and that it is probably a threshold phenomenon, simi- 
lar to a second-order phase transition. Similar results 
were earlier obtained in the field experiment by Banner, 
Babanin and Young 12 1. 



DYNAMICAL MODEL 

In this part of our experiment, the surface of the liquid 
is described by two functions of the horizontal variables 
x,y and the time t: the surface elevation rj(x,y,t) and 
the velocity potential on the surface ip{x,y,t). In our 
approximation, they satisfy the following equations Q: 
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Here k is the linear integral operator k = v 
corresponds to the inverse Fourier transform. 

Equations |T]) are nowadays widely used in numerical 
experiments and are solved by different versions of the 
spectral code [H, Q G3, ES El 0, El, M, 0, 0, [24 1 . 
In the present experiment, we solved the equations in 
the real space domain 2tt x 2tt using the finest currently 
possible rectangular grid 512 x 4096, putting g = 1. The 
dissipative terms F 1-1 ^?^] and F -1 ^^^] are taken in 
the form of pseudo- viscous high frequency damping. We 
put 



7fc 
kd, 



0,k< k d , 

-j{k - k d ) 2 ,k > k d 
1024,7 = 5.65 x 10~ 3 . 
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In accordance with recent results 25[, the dissipation 
term should be included in both equations. 

The distribution of the wave action is described by the 
function n(k,t) = |ar(t)| 2 , where 



(3) 
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are complex normal variables. Here u>k = \fgk. 

As the initial condition, we used a Gaussian-shaped 
distribution in the Fourier space: 
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(4) 

The initial phases of all harmonics were random. The 
average steepness of this initial condition, defined as \x = 
v/2(|Vr/| 2 ), was fi ~ 0.176. 

The period of the most intensive wave was To = 
27r/v / 300 = 0.362. Calculations continued until t = 
3378T . We observed an angular spreading of the ini- 
tial spectral distribution together with a downshift of 
the spectral peak. Level-lines of the initial and the fi- 
nal spectra are presented on Figs. [TJ [U We observed 
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FIG. 1: Initial spectrum |a^| 2 . t = 0. 
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FIG. 2: Final spectrum |a^| 2 . t ~ 2809T . 



the following indications of wave- turbulent behavior: 
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1. The statistics of energy-capacity spectral modes is 
close to the Rayleigh distribution. We observed 
the presence of a few very intensive harmonics (so- 
called oligarchs, cf. which did not obey the 
Rayleigh statistics, but their contribution to the 
total balance of the wave action is small (no more 
than 5%). This means that we almost overcame 
negative effects caused by the finite size of our sys- 
tem (see 21, 22, HH), and that our grid is fine 
enough. 

2. We observed the formation of the Zakharov- 
Filoncnko spectral tail in the energy spectrum 2 
(see Fig. El). 




FIG. 3: The time-averaged spectrum of a point on the surface 
n u = (|j7u| 2 ) in the double logarithmic scale. The tail of the 
distribution fits to the Zakharov-Filonenko spectrum uj~ 4 . 

At the same time, we observed a manifestation of strong- 
turbulent effects. They are manifested by the formation 
of "fat tails" on the PDF for surface elevations and espe- 
cially for its gradients (see FigSJ) The presence of these 
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FIG. 4: PDF for (Vn) y at the moment of the maximum sur- 
face roughness, t ~ 14To. 

tails indicates that the surface has a tendency to become 
rough and to produce white-capping. In our model, wave- 
breaking is arrested by the strong pseudo- viscosity. 



STATISTICAL EXPERIMENTS 



In the second experiment, we solved the Hasselmann 



kinetic equation for n? = 
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Here jk is the pseudo-viscosity and Sdiss is the phe- 
nomenological dissipation term modelling the white- 
capping process. 

Eq. §5§ was solved on the grid 71 x 36 in polar coor- 
dinates on the frequency-angle plane by the Resio- Tracy 
code , improved in [I, @| ■ We first performed the ex- 
periment with Sdiss = 0. We observed good qualitative 
coincidence with the dynamical experiment. We observed 
a downshift of the spectral peak, angular spreading and 
the formation of o-> -4 spectral tails. But the quantitative 
agreement of the experiments was not good: it was clear 
that the inclusion of some phcnomenological dissipation 
is necessary. 

We examined the standard from of Sdiss used in in 
the industrial operational models of wave forecasting — 
WAM Cycle 3 and WAM Cycle 4 (hereafter WAM3 and 
WAM4) [26|: 



Sdi. 



k 
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where k and to are the wave number and the frequency, 
tilde denotes the mean value; Cd s , $ and p are tun- 
able coefficients; S — k\/~H is the overall steepness; 
Spm = (3.02 x 1CT 3 ) 1 / 2 is the value of § for the Pierson- 
Moscowitz spectrum (note that the characteristic steep- 
ness is /i ~ V%S). It is worth noting that according to 
[I2I ] , the theoretical value of the steepness for the Pierson- 
Moscovitz spectrum is Spm — (4.57 x 10 -3 ) 12 , which 
gives us fi ~ 0.095. 

The values of tunable coefficients in the WAM3 case 
are: 

Cds - 2.36 x 10~ 5 , 5 = 0, p = 4 (7) 

and in the WAM A case are: 

C ds = 4.10 x 10~ 5 , (5 = 0.5, p = A (8) 

The evolution of the total wave action is presented on 
Fig. [5] One can see that in the long run, the models 
WAM3 and WAM4 overestimate white-capping dissipa- 
tion. To achieve better agreement of both experiments, 
we used the following form of the dissipative term: 

C ds = 1-00 x 10~ 6 , 5 = 0, p=12 (9) 
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The total wave action curve corresponding to this new 
dissipation term is shown on Fig. [5] by the thick solid line 
and displays excellent correspondence with the dynami- 
cal model. 
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FIG. 5: Total wave action as a function of time. The solid line 
corresponds to the dynamical equations, the dashed-dotted 
line — to the kinetic equation with artificial viscosity, the 
dashed line — to the kinetic equation with the WAM3 damp- 
ing term, the dotted line — to the kinetic equation with the 
WAM4 damping term, and the thick solid line — to the ki- 
netic equation with the new damping term, 
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CONCLUSION 

Our experiments can be interpreted as a confirmation 
of the theory of weak turbulence. However, even at mod- 
erate values of the parameter of the nonlinearity //, the 
strongly nonlinear effects of white-capping are essential. 
They manifest themselves as fat tails of the PDF and 
lead to additional dissipation of wave energy. This dis- 
sipation demonstrates a very strong dependence on the 
steepness. At steepness fi = 0.176 they dominate, at 
steepness fx = 0.09 they are negligibly small. The re- 
sults of our experiments are in good qualitative agree- 
ment with the field experiment of Banner, Babanin and 
Young [l2] , who found that wave- breaking is a threshold 
effect, similar to a second-order phase transition. 

We stress that the dependence (|9]) is much sharper than 
it is usually stated. So far, the sharpest dependence p = 5 
was given by Donelan [28| . We can guess that the real de- 
pendence of Sdiss on \x is even stronger, and that the on- 
set of the wave breaking is a threshold-type phenomenon 
like a second-order phase transition. 
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